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In this paper, we consider three-dimensional dynamical systems, as for example the Lorenz 



model. For these systems, we introduce a method for obtaining families of two-dimensional 
surfaces such that trajectories cross each surface of the family in the same direction. For 
obtaining these surfaces, we are guided by the integrals of motion that exist for particular 
values of the parameters of the system. Nonetheless families of surfaces are obtained for 
arbitrary values of these parameters. Only a bounded region of the phase space is not filled 
by these surfaces. The global attractor of the system must be contained in this region. In 



this way, we obtain information on the shape and location of the global attractor. These 
results are more restrictive than similar bounds that have been recently found by the method 
of Lyapunov functions. 
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The Lorenz equations ([]]) are one of the classic models of nonlinear dynamics and 
chaos. These equations were originally derived in a modal truncation of the Boussinesq 
equations for thermal convection. They read as follows : 

x = a(y — x) 

y = rx — y — xz (1) 

z = xy — bz 

with cr,b,r > 0. There a corresponds to the Prandtl number, b is a geometric parameter 
and r is the Rayleigh number in units of the critical Rayleigh number. 

These equations describe a dissipative dynamical system for all values of r, a and 
b because the divergence of the flow field is always negative. Hence 3-dimensional 
volumes in the phase space contract to zero at a uniform exponential rate and the 
system's attractor is necessarily of dimension less than three. This model has become 
a classic in the area of nonlinear dynamics. Its importance is not that it quantita- 
tively describes the hydrodynamics motion, but rather that it illustrates how a simple 
model can produce very rich and varied form of dynamics, depending on the value of 
a parameter in the equations ||. 

In this paper, we are interested in the approximated location in the phase space 
of the global attractor of the system, which contains all dynamics evolving from all 
initials conditions. The global attractor is the set of points in phase space that can be 
arrived at from some initial condition at an arbitrary long time in the past. The two 
fundamental properties of global attractors are || : 

• it is invariant under the evolution. 

• the distance of any solution from it vanishes as t — > +oo. 

The last property is simply interpreted thus : if the solution starts initially outside 
the global attractor, then it is attracted into it as t — > +oo and once inside it cannot 
escape. If it starts inside then it stays inside. 

The global attractor contains all the asymptotic motion for the dynamical system. 
It is common to talk of multiple attractors for a dynamical system, and each of them 
may in its own right be considered as the attractor for initial conditions within its 
own bassin of attraction. The notion of global attractor corresponds to the union of 
all possible such dynamically invariant attracting sets. In particular, it contains all 
possible structures such as fixed points, limit cycles etc... 



The global attractor is contained in an absorbing ball in phase space, and we want 
to obtain analytic estimates about its geometric shape. Moreover, this enables us to 
find good estimates of its Lyapunov dimension. Estimates which give the shape of the 
attractor are important as they lead to a good upper bound on the dimension of the 
Lorenz attractor ||. 

Until now, approximate locations of the Lorenz's attractor in the phase space have 
been obtained by the method of Lyapunov functions [0, |, || ||, || . Very recently, thanks 
to this method, it has been shown that the global attractor of the Lorenz equations 
is contained in a volume bounded by a sphere, a cylinder, the volume between two 
parabolic sheets, an ellipsoide and a cone . 

In this paper, we apply a different method for obtaining analytic estimates for the 
location and shape of the Lorenz attractor. The method is based on the determination 
of families of 2-dimensional surfaces that are crossed by the trajectories of the system 
only in one direction. In the region filled by these surfaces, the dynamical behaviour 
is very simple. The asymptotic complex behaviour must be contained in the region of 
the phase space that is not occupied by these surfaces. 

For finding these families of surfaces, we will be guided by the time-dependent inte- 
grals of motion that exist for special values of the parameters of the system. Integrals 
of motion for the Lorenz system have been extensively studied in p], [|, [7| . The known 
integrals of motion are : 

a I(x, y, z, t) = (x 2 — 2az)e 2at with b = 2a and a and r arbitrary. 

b I(x, y, z, t) = (y 2 + z 2 )e 2t with b = 1, r = and a arbitrary. 

c I{x, y, z, t) = (—r 2 x 2 + y 2 + z 2 )e 2t with b = 1, a = 1 and r arbitrary. 

d I(x, y, z, t) = (^^x 2 + ay 2 - (4a - 2)xy - ±x A + x 2 z) e 4<rf 
with b = 6a — 2, r = 2a — 1 and a arbitrary. 

e I(x, y, z, t) = (— rx 2 — y 2 + 2xy + \x 4 — x 2 z + 4(r — 1)2:) e 4t 
with b = 4, a = 1 and r arbitrary. 

For each of these integrals we have % = 0. Let us consider case a and let us define the 
family of surfaces x 2 — 2az = k, where k is an arbitrary constant. The scalar product 
between the normal vector N to this surface at a given point and the tangent vector 



T to the trajectory of the Lorenz system that goes through this point is given by : 

N • T = (2xi — 2ak) ■ (a(y — x)i + (rx — y — xz)j + [xy — bz)k\ = —bk (2) 

Therefore, for a given surface (i.e. for a given value of k) this scalar product has the 
same sign for all the points of the surface. Each surface of the family is crossed in the 
same direction by the flow associated to the system. This direction depends of the sign 
of the constant k. Hence, for the case b = 2a, the 3-d phase space of the Lorenz system 
is filled by two families of surfaces, the families associated to positive and negative 
values of k. The scalar product N ■ T is positive (resp. negative) for negative (resp. 
positive) values of k. It is clear that the surface corresponding to k=0 plays a very 
special role. All the trajectories of the system are attracted by this surface. On this 
surface, the scalar product N ■ T is zero. This surface is an invariant manifold of the 
system, as can be seen in fig. [3]. 
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Figure 1: The family of surfaces x 2 — 2oz = k for the case b = 2a. The bolded surface 
corresponds to k = and is the attracting set of the system. Some trajectories of the 
system are shown. 



It is clear that the existence of these families of surfaces gives a lot of information 
about the dynamics of the system. The behaviour of trajectories is extremely simple 
in all the phase space with the exception of the invariant surface x 2 — 2a z = 0. This 
surface contains the global attractor of the system for the case b = 2a. Here, the global 
attractor is contained in a two-dimensional surface, as for the five cases a, b, c, d and 
e, that is when an integral of motion exists. The family of surfaces derivated above 
enables us to characterize in a simple way this global attractor. The determination of 



this family of surfaces follows immediatly from the existence of the integral of motion 
a, when b = 2a. Now, the natural question is : if b ^ 2a, is it possible to find similar 
families of surfaces that the flow crosses in the same direction at each point of the 
surface ? - in the following, we will call this type of surface semipermeable - In this 
case, we do not have at our disposal an integral of motion, and these surfaces can not 
fill the phase space because, in the general case, the global attractor is not contained 
in a two-dimensional set. 

In order to find semipermeable surfaces in the general case (when integrals of motion 
do not exist), we will procede as follows : 

For the case of the integral of motion a, we first propose, when b ^ 2a, a surface of the 
same mathematical form as the integral of motion a, but with arbitrary coefficients : 

S = a\x 2 + a2Z + a 3 = (3) 

The scalar product N ■ T is now : 

N ■ T = 2a\xx + a 2 z 

= {2aia + a 2 )xy — 2a\OX 2 — a 2 bz (4) 

If we calculate this scalar product on the surface S, for the general case b ^ 2a, we 
obtain : 

N ■ T/S = (2oicr + a 2 )xy + (ba\ — 2<ra 1 )x 2 + ba 3 (5) 

where we have replaced in (f|) —02-2 by 03 + a\x 2 . We now have an expression that 
depends only on two variables : x and y. The problem of determining the coefficients 
Oi, a 2 and 03 in order for this expression to have the same sign for arbitrary values of 
x and y is considerably simpler than the analogous problem in the three variables x,y 
and z that must be solved in the method of the Lyapunov functionals. 

To keep the same sign in (|5[) for arbitrary values of x and y, we must take a 2 = 
—2aa\. Then we have : 

N ■ f/S = 01 (b - 2a)x 2 + ba 3 (6) 

As a\ must be nonzero, we can take ai = 1 without loss of generality. We now have 
two different cases : 

i) b > 2a, we must take 03 > in order to have a family of semipermeable surfaces. 
We show this family, as well as some trajectories of the system, in fig. S 



ii) b < 2a, we must take a^ < in order to have a family of semipermeable surfaces. 
We show this family, as well as some trajectories of the system, in fig. ||]. 




Figure 2: The family of semipermeable surfaces (|3]) for the case b > 2a. The bolded 
surface is the last surface of the family. Some trajectories of the system are also shown. 
The critical points C+ and C- are below these surfaces. 

As we can see from the figures above, in the region filled by the surfaces the dynamic 
of the system is very simple. The complex behaviour can only occur in the region of 
phase space that is not occupied by these surfaces. In case ii), the global attractor 
of the system must be located in the region z > 2a x 2 . For the case i), because of 
the presence of the semipermeable surfaces, the flow cannot enter the z > 2ax 2 region 
upward and hence the homoclinic trajectory cannot exist. 

Therefore, motivated by the existence of the first integral a, valid in the case b = 2a, 
we have found a family of semipermeable surfaces for arbitrary values of the parameters 
of the system. 

As we shall see, new families of semipermeable surfaces can be found by using the 
other integrals of motion. 



From the case b, we deduce that the surfaces y 1 



k 2 are semipermeable for 



arbitrary values of k, when b = 1 and r = 0. Guided by this result, we propose in the 
general case a family of surfaces of the form : 



S = ai(y - ci) 2 + a 2 (z - c 2 ) 2 -1 = 



(7) 




Figure 3: The family of semipermeable surfaces^) for the case b < 2a. The bolded 
surface is the last surface of the family. Some trajectories of the system are also shown. 
The chaotic attractor is above these surfaces. 

The scalar product N • T is given by : 

N ■ T = 2((a 2 — di)xyz — a 2 bz 2 — a±y 2 + (a^ — a 2 C2)xy + a\C\XZ 

+a 2 c 2 bz + aiCiy — a\C\rx) (8) 

The evaluation of N ■ T on the surface (|7]) becomes far simpler if we take a 2 = a\, 
c\ = and c 2 = r. After this, (0) and (§) respectively become : 



S = ai(y 2 + (z-r) 







N • f = -la^z 1 + y 2 - rbz) 
The scalar product iV • T calculated on the surface S is given by : 



N ■f/S = a l ({\- b)z 2 + r(b - 2)z + r 2 - 



a i 



(9) 
(10) 

(11) 



Note that N-T/S is a function of only one variable, as it is the case for expression (H) 
for the semipermeable parabolas. In this case, the surface S is not infinite in the y and z 
directions. In particular, the coordinate z varies in the interval : r — -== < z < r + -^=. 
In consequence, the quadratic polynomial (|Tl|) must have the same sign only in this 
interval and not for arbitrary values of z. This condition determines the possible values 
of <2i, that can be found by applying the Sturm's theorem. The results are as follows : 



b < 2 , — > r 2 

1 r 2 b 2 

b> 2 , — > — 

' a x ~ 4(6 - I 



(12) 



Hence, for arbitrary values of the parameters of the system, we have found a family of 
semipermeable infinite cylinders. The radius of these cylinders varies between +00 and 
the minimal values given in fll2|) . The behaviour of some orbits with respect to this 
family of surfaces is shown in fig. f|. The global attractor is contained in the region 
not occupied by these surfaces. 




Figure 4: Chaotic attractor stuck inside semipermeable cylinders. The bolded circle 
corresponds to the smallest cylinder. 

From the integral of motion c, we deduce the existence of the family of semiper- 
meable surfaces z 2 + y 2 — rx 2 — k = 0, with k arbitrary and b — 1, a — 1. 
Guided by this result, we propose in the general case the family of surfaces : 

S = a\X 2 + a 2 y 2 + a 3 z 2 — a 4 = , with 01 a 2 a 3 < (13) 

The scalar product N • T is given by : 

N ■ T = 2 ((03 — a 2 )xyz — a 3 bz 2 — a 2 y 2 — a\ox 2 + (oicr + a 2 r)xy) (14) 

In order to have the same sign for N -T on the surface S, we take a 3 = a 2 = 1. Then 
we have : 

S = aix 2 + y 2 + z 2 — a 4 = , with a\ < (15) 

N'T = 2 (-(y 2 + bz 2 ) - a x ax 2 + [a x o + r)xy) (16) 

This expression, calculated on the surface (|T5|), gives : 

N ■ f/S = 2 (01(6 - o)x 2 + (aicr + r)xy + (b - l)y 2 - ba A ) (17) 
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On the surface flT5|) the variables x and y vary in such a way that the following inequality 
must be satisfied : a\x 2 + y 2 < a^. Therefore (|17|) must have the same sign for all 
values of x and y that satisfy this inequality. The solution of this algebraic problem is 
not simple. Hence we do not give here the technical details of the calculations. 
This condition determines the possible values of the coefficient a\ : 

-2or - (a - l) 2 - y/(a-iy + 4ar(a-l) 2 



-2ar 



<J 



2a 2 
-l) 2 



< a x < 
v/(a-l) 4 + 4ar(a-l) 2 



2a 2 



or 



2(6- 


- l)(b-a)- ar -2^/{b- l) 2 (fo- 


- a) 2 — ar{b - 


!)( & -^/„ / 




O 2 








2(b - l)(b - o) - or + 2^{b~- 


-\) 2 {b-a) 2 


-ar(b- l){b-a) 



<y- 



(19) 



The parameter a 4 is arbitrary, and the condition a\ < restricts the possible values of 
the parameters b and a, which are given in fig. |5|. The canonical values r = 28, a = 10 




Figure 5: Parameters a and b for which surfaces (|T5| ) are semipermeable. 

and b = | lie in the region I. For parameters a and b in region I of fig. |5|, there are two 
zones in phase space that are filled with surfaces of the family. The behaviour of some 
trajectories of the system, in relation to the semipermeable surfaces is shown in fig. [| 

These results are more restrictive (they give a more precise information about 
the location of the global attractor) than similar results obtained recently in by 
employing the method of Lyapunov functions. Therefore, our method can locate more 
accurately the global attractor of the system in phase space than before. 

In region II of fig. |^, the scalar product iV • T/S has opposite sign with respect to 
region I and, in phase space, there are two zones that are not filled by the surfaces of the 
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Figure 6: Lorenz attractor squeezed between semipermeable hyperboloids ([T5|) . The 
bolded lines are the last repelling cones. The thick dashed line is the former bound 
obtained by Doering & al. ||. The parameters a and b are in region I. 

family. These two zones are not connected between them and the critical points C- and 
C+ are contained in each of the different regions. The behaviour of some trajectories 
of the flow with respect to the semipermeable surfaces and the position of the critical 
points C- and C+ are shown in the fig. [7]. If one trajectory enters one of the two free 
regions in the phase space, it cannot exit from it and hence cannot pass in the other 
one. This restriction on the behaviour of the orbits prevents the possibility of a chaotic 
behaviour. The trajectories that evolve around one of the critical points cannot go to 
the other free region for evolving around the other critical point (we make reference 
to the critical points C- and C+ only). It is clear that the homoclinic bifurcation that 
precedes the birth of the chaotic behaviour cannot occur in region II of the parameter 
space. 

There is still another result for region II of the parameter space : since the flow, once 
it has entered one of the two free regions of the phase space, cannot escape from it, it 
can only go to the critical point C- or C+ lying in this region. So each one of the two 
free zones in the phase space is a part of the attraction's basin of C- or C+. 

By applying the same method and guided by the form of the integral of motion 
d, we have found a new family of semipermeable surfaces. We do not give here the 
technical details of the calculations. They are a little more complicated than the 
calculations involved in the previous cases. The results are as follows. The family of 
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Figure 7: Critical points C- and C+ separated by semipermeable surfaces (O). The 
left (resp. right) zone is a part of the bassin of attraction of C- (resp. C+). The 
parameters a and b are in region II. 

surfaces : S = c\ + c^x 2 — -^x 4 + (2a — b)xy + cry 2 + x 2 z = are semipermeable in the 
two following cases : 

i) 6 < 6er — 2, r < 2er — 1, Ci > and c<i G I C2 , where I C2 is the interval defined by 
the two real roots of the quadratic polynomial in c<i : Aa 2 c\ + 4a (b + 2er — bo + 
2ra - 6a 2 )c 2 + b 2 - Aba - 6b 2 a - Abra - Ab 2 ra + Ao 2 + 2Abo 2 + % 2 a 2 + 8ra 2 + 
20bra 2 + Ar 2 a 2 — 24<r 3 — 36&<r 3 — 2Ara 3 + 36a 4 . In this case, the sign of the scalar 
product N ■ T/S is positive. The canonical values of the Lorenz's parameters do 
not satisfy the above two inequalities between r, b and a. 

ii) b > 6a — 2, r > 2a — 1, c\ = and c 2 G I C2 . 

In this case, the sign of iV ■ T/S is negative. This family of surfaces divides the 
phase space in three regions. Only one of them is filled by the surfaces of the 
family. The two free regions are disconnected and the critical points C+ and Cl- 
are located in different regions. Here, as in one of the cases analysed above, the 
homoclinic bifurcation, and hence the chaotic behaviour, is not possible. The 
critical points C+ and C- are stable and each of the two free regions is part of 
the basin of attraction of each critical point. 

Finally, guided by the form of the first integral e, we have found another family of 
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semipermeable surfaces. Let us consider the family of surfaces : 



S = Cl + 



2 - 36 + b 2 - c 2 + 2a - ba - 2ra 



-x 



2a Aa 

These surfaces are semipermeable in four different cases : 



—x 4 + (2 - b)xy + ay 2 + (c 2 + x 2 )z = 



(20) 



a b> 2{a + 1), a < 1, Cic 2 < 0, c x < 



(c2+k3)(c2k2+b 2 k:j) , 
64(cr-l)cr(o-+l) ' 



here the flow is crosssing the surfaces downward. 

13 b>2(a + 1), a < 1, Cl c 2 < 0, c 2 > ^fg^, Cl < 2-36+^-^+2.-^-^ . 
here the flow is crossing the surfaces downward too. 

7 6 < 2(0- + 1), cr > 1, c lC2 > 0, Cl > ^ggg^ ; 

here the flow is crossing the surfaces upward. 



5 b < 2(cr + 1), a > 1, cic 2 > 0, (6 - 4a + 4)c 2 < -bk 3 

C l > 2^ ' 

here the flow is crossing the surfaces upward too. 



The quantities fc 2 and /C3 are given by : k 2 = b 2 + 8b(a — 1) + 16(1 — a 2 ) and 
k 3 = 2b - b 2 - Aa + 26a + Ara. 




Figure 8: Critical points separated by semipermeable surfaces of type (BO) in case a. 
The far left (resp. right) zone is a part of the bassin of attraction of C- (resp. C+). 



If c 2 > 0, each surface of the family is connected; if c 2 < 0, all the surfaces are 
disconnected : they are divided in three parts. 
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Figure 9: Lorenz attractor enclosed by semipermeable surfaces of type fl2"UD in case 7. 
The form of these surfaces for x 2 < — c 2 is a sink, what we cannot see in the figure, 
which is a projection. The bolded curve is given by the equality in expression ([H|). 
The thick dashed curved is the last parabola. 

In case a, if we take r > ■^(b — 2)(b — 2cr), c 2 G [— b 2 -^-; — k 3 ] < and then c\ > 0, we 
then have disconnected semipermeable surfaces and the critical points are under the 
surfaces and separated by them. This is another configuration where we know a part 
of the basin of attraction of each critical point and where the homoclinic trajectory 
cannot exists (see fig. |8|). 

The case 7 gives us information about the space extension of the chaotic attractor 
(r = 28, a = 10, b = §). In this case, when c 2 < and c\ > 0, the surfaces are 
disconnected and the flow crosses them upward (see fig. [5J). The uppermost surface 
(for c\ = and c 2 = — ks) is an additional bound for the Lorenz attractor. Hence it 
lies entirely in the zone of the phase space where : 



z > 



2-36+6 2 - C2 +2 g -6 g -2rcr a ,2 _ J_ x 4 + q _ ^ + ay l 

x 2 — k? 



(21) 



(for Lorenz's canonical values ^3 = 1131.56). 

We have also found several semipermeable families of ellipsoids. In fact, we have 
generalised results given in || [|, ||. Surfaces like 



S =^^l X 2 + y 2 + ( Z - C3 ) 2 = R 



a 



(22) 



are semipermeable for the following cases : 
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R>c; ° (>cp 
4(b-l) 



R >c 



Figure 10: Parameters for which surfaces fl22|) with C3 > r are semipermeable 



C3 > r for arbitrary a, b, r and for values of R as in fig. [nj. The interest of 
having one free parameter (here C3) in addition to R is that we may lower the 
ellipsoids in phase space with C3 and so restrict more tightly the region in which 
the chaotic attractor lies (with considering the envelope of all the smallest (with 
R minimum) ellipsoids when C3 varies in ]r; + oo[). These results contain known 
results about ellipsoids and several new ones. 

if C3 = r then S is the cylinder which we have studied above. 

if C3 < r then S is an hyperboloid of revolution. The revolution axis is {z = 
C3,y = 0}. This surface is semipermeable for b < a and a > 1 (a case which 



includes Lorenz's canonical values) and R < 



b 2 



Ao(b-o) 



C3 2 For this case the scalar 



product is positive V {x, y} G S. The flow crosses the last surface 



c^—r 2 



x* +y* + 



{* - c 3 



b 2 



4cr(b-a) 



C3 2 outwards. This new information sharpens the bounding 



frontiers of the attractor as we can see in fig. 1 1 . 



Sparrow || has conjectured that all trajectories of the Lorenz system eventually 
enter and remain in the region z > for all parameters values r, a and b (note that 
the plane z = is not semipermeable). Sparrow proved this conjecture for the case 
b < a + 1 by using the method of Lyapunov functions. 

2 

The existence of semipermeable parabolas z = |- + ao for ao < and b < 2a also 
proves the conjecture, but for different values of parameters a and b. Indeed, in this 
case the flow is crossing all the parabolas upward. It doesn't mean that the z = 
plane itself is semipermeable, but the flow has to end up with crossing upward the last 

2 

parabola z — f-. The parabolas are in fact driving the flow to the phase space zone 
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Figure 11: Lorenz's attractor and a semipermeable surface of the family fl2"2"|). When we 
consider all the surfaces of the family, stronger bounds on the location of the attractor 
are obtained. 



2 

where z > §- > 0. 



The cylinders family : y 2 + (z — r) 2 = R with R> r 2 and b < 2 also drives the flow 
inside the smallest cylinder (R = r 2 ) which lies in the z > phase region. 

The surfaces (p0|), in case 7, tell us also that the flow eventually crosses the up- 
permost surface (given by the equality in (|2l|)) upward. This surface lies entirely in 
the z > half space (for x 2 < — C2). The results of this work expand the region of 
the parameter space for which the flow eventually enters the zone of phase space with 



z > 0. This region is shown in fig. 12. 



In conclusion, inspired by the integrals of motion that exist for particular values of 
the parameters b, o and r, we were able to find several families of surfaces, all crossed 
in the same direction by the flow associated to the system. 

From these results, we have deduced a rich quantity of information about the geo- 
metrical location of the global attractor of the system. This information is more restric- 
tive than similar results that had been found by the method of Lyapunov functions. 
When compared to the Lyapunov technique, we see that the fundamental advantage 
of this new method is that one now has to study functions with one less variable. 

Moreover, we have obtained information about the spread of the basin of attraction 
of the critical points C- and C+ when they are stable. We have also determined regions 
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Figure 12: Range of parameters a and b for which the flow eventually enters the phase 
space zone z > (Vr) 

of the parameter space where the chaotic behaviour is not possible. 

It is clear that the method used in this paper can be applied to other 3-d dissipative 
dynamical systems that the Lorenz one. We have choosen the latter owing to the great 
importance that this system has played in the study of chaotic dynamics. 
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